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Lecture notes 



Lecture 8 



In this lecture, the electric polarisation density of the medium is finally inserted into Maxwell's 
equations, and the wave propagation properties of electromagnetic waves in nonlinear optical media 
is for the first time in this course analysed. As an example of wave propagation in nonlinear optical 
media, the optical Kerr effect is analysed for infinite plane continuous waves. 

The outline for this lecture is: 

• Maxwells equations (general electromagnetic wave propagation) 

• Time dependent processes (envelopes slowly varying in space and time) 

• Time independent processes (envelopes slowly varying in space but constant in time) 

• Examples (optical Kerr-effect xj^^^^(— w; w, w, — w)) 

Wave propagation in nonlinear media 

Maxwell's equations 

The propagation of electromagnetic waves are, from a first principles approach, governed by the 
Maxwell's equations (here listed in their real- valued form in SI units), 

dB(r t) 

V X E(r, t) = (Faraday's law) 

dT>(r t) 

V X H(r, t) = J(r, t) H (Ampere's law) 

V-D(r,t) =p(r,<), 
V-B(r,t) = 0, 

where p{r,t) is the density of free charges, and 3(v,t) the corresponding current density of free 
charges. 

Constitutive relations 

The constitutive relations are in SI units formulated as 

D(r,i) =£oE(r,i)+P(r,i), 
B(r,t) = /^o[H(r,t) + M(r,t)], 

where P(r, f) = P[E(r, t), B(r, t)] is the macroscopic polarization density (electric dipole moment 
per unit volume), and M(r. f) = M[E(r, t), B(r, f)] the magnetization (magnetic dipole moment 
per unit volume) of the medium. 

Here E(r, t) and B(r,f) are considered as the fundamental macroscopic electric and magnetic 
field quantities; D(r,f) and ll{r,t) are the corresponding derived fields associated with the state 
of matter, connected to E(r, <) and B(r, t) through the electric polarization density P{r,t) and 
magnetization (magnetic polarization density) M(r, t) through the basic constitutive relations. In 
fact, the constitutive equations above form the very definitions^ of the electric polarization density 
and magnetization. 



^ J. D. Jackson, Classical Electrodynamics, 2nd cd. (Wiley, New York, 1975); J. A. Stratton, 
Electromagnetic Theory (McGraw-Hill, New York, 1941). 
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Two frequent assumptions in nonlinear optics 

• No free charges present, 

p(r,t) = 0, J(r,i)=0. 

(Any relaxation processes etc. are included in imaginary parts of the terms of the electric 
susceptibility.) 

• No magnetization of the medium, 

M.{T,t) = 0. 



The wave equation 

By taking the cross product with the nabla operator and Faraday's law, one obtains 

d 

V X V X E(r,t) = -^V X B(r,t) 
d 

= -/«Og^VxH(r,t) 
d dT>{T,t) 

/ d'^{r,t) ^'P{r,t) ^ 

Since now hqSq = 1/c^ in SI units, with c being the speed of light in vacuum, one hence obtains 
the basic wave equation, taken in time domain, as 

VxVxE(r,0 + ^-^ = -Mo^;^, (1) 

where, as in the previous lectures of this course, the polarization density can be written in terms 
of the perturbation series as 

P(r,t) = ^pW(r,i)=£o J xW(-u;;u;)iJ„(r,u;)exp(-za;t)du; + ^p('=)(r,t) 

=P(i>(r,t) =p(NL)(r,t) 

In the left hand side of Eq. (1), we find the part of the homogeneous wave equation for propagation 
of electromagnetic waves in vacuum, while the right hand side described the modifications to the 
vacuum propagation due to the interaction between light and matter. In this respect, it is now 
clear that the electric polarisation efiiectively acts as a source term in the mathematical description 
of electromagnetic wave propagation, making the otherwise homogeneous vacuum problem an 
inhomogeneous problem (though with known source terms). 

It should be noticed that whenever the polarization density is calculated from the Bloch equa- 
tions (formulated later on, in lecture 10 of this course), instead of by means of a perturbation 
series as above, the Maxwell equations and the wave equation (1) above are denoted Maxwell- 
Bloch equations. In some sense, we can therefore see the choice of method for the calculation of 
the polarization density as a switch point not only for using the susceptibility formalism or not for 
the description of interaction between light and matter, but also for the form of the wave propa- 
gation problem in nonlinear media, which mathematically significantly differ between the "pure" 
Maxwell's equations with susceptibilities and the Maxwell-Bloch equations. 
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The wave equation in frequency domain (optional) 

Frequently in this course, we have rather been studying tlie {^kx;tric fields and polarisation densities 
in frequency domain, since many phenomena in optics are properly and conveniently described as 
static (in which case the frequency dependence is simply reduced to the interaction between discrete 
frequencies in the spectrum). By using the Fourier integral identity^ 



/oo 
E^{uj) exp(-iu;i) dio = ^-^[E„]{t), 
-oo 



with inverse relation 



1 

Ea{^^) = J ^oc{t) exp(iwr) dr = ^[Ea\{uj), 
we obtain the wave equation (1) as 

^2 

V X V X E(r,a;) - -^E(r,a;) = iiou?'P{y,uj). 
Quasimonochromatic light - Time dependent problems 

By inserting the perturbation series for the electric polarisation density into the general wave 
equation (1), which apply to arbitrary electric field distributions and field intensities of the light, 
one obtains the equation 

1 92 a2p(NL)/j. ^\ 

V X V X E(r,t) + ^— / e^e^a{'^)Ea{r,uj)eyip{-iujt)dbJ = -iiQ — (2) 



where 



(denote this integral as / for later use) 



is a parameter commonly denoted as the relative electrical permittivity.^ This wave equation is 
identical to Eq. (7.14) in Butcher and Cotter's book. (Notice though the printing error in Butcher 
and Cotter's Eq. (7.14), where the first ^0 should be replaced by 1/c^.) 

The second term of the left hand side of Eq. (2) gives all first order optical contributions to 
the wave propagation, as well as all linear optical dispersion effects. This terms deserves some 
extra attention, and we will now proceed with deriving the effect of the frequency dependence of 
the relative permittivity upon the wave equation. First of all, we notice that since £'a(r, — w) = 
E'* (r, Lo), which simply is a consequence of the choice of complex Fourier transform of a real valued 
field, the reality condition of Eq. (2) requires that 



^ From the inverse Fourier integral identity, it follows that the Fourier transform of a derivative 
of a function f{t) is 

^[/'W]H = -^'^^[/(i)]M => ^[/"(t)](c^) = -c.2 5[/(i)]H- 



^ Notice that for isotropic media, ■x^^a{—uj\oj) = Xxx (— w)(5^a) which leads to the simplified 
form 

Gn£^a{i^)Ea{r,Lo) = £(a;)E(r,a;). 

We will here, however, continue with the general form, in order not to loose generality in discussion 
that is to follow. 
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It should be emphasized that this property of the relative electrical permittivity merely is a conve- 
nient mathematical construction, since we in regular physical terms only consider positive angular 
frequencies as argument for the refractive index, etc. 

For quasimonochromatic light, the electric field and polarisation density are taken as 



E(r,i)= Re[E^^(r,i)exp(-iw^t)], 
P(r,i)= Re[P„^(r,i)exp(-ia;^t)], 



a;„>0 



where E^^^^ (r, t) and P^^^ (r, t) arc slowly varying envelopes of the fields. In the frequency domain, 
the quasimonochromatic fields are expressed as 



E(r, ^) = \Y. (r, ^-^.)+ (r, - a;<.)], 



2 

where the envelopes have some limited extent around the carrier frequencies at ztcjcr- Notice that 
the fields taken in the frequency domain are expressed entirely in terms of their respective temporal 
envelope, that is to say, without the exponential functions that appear in their counterparts in time 
domain. 

For simplicity considering a medium that in the linear optical domain is isotropic, with the 
relative electrical permittivity 

£^ta{'^) = e(w)5^a = nl{ijj)6u,a, 

where no(ijj) is the first order contribution to the refractive index of the medium, this leads to the 
middle term of the wave equation (1) in the form 



/oo 
e^e/iQ(w)S„(r,u;) exp(-iwf) dhj 
-OO 

/ ^ ^ V [E^„ (r, a; - w^) + E* ^ (r, -u - to^)] exp{-iLjt) dio 



quasimonochromatic form of E(r,aj) 

{denote uj"^ e{Lo) / (? = uj'^n]^{Lo)/c^ = fc^(w)} 



oE / fc^H[E^<,(r,a;-w<,) + E*^(r, -w-w<,)]exp(-iwt)dw. 

1 r°° 

- Y^ / fc^(w)E;j_^(r, w — Wct) exp(— iwi) + c. c. 



If now the field envelopes decay to zero rapidly enough in the vicinity of the carrier frequencies 
(as wc! would expect for quasimonochromatic; light, with a strong spectral confinement around the 
carrier frequency of the light), then we may expect that a good approximation is to make a Taylor 
expansion of fc^(w), in the neighbourhood of respective carrier frequency of the light, as 



dk 
dL 

dk 

duj 



^ kfj -\- 2/C(7 - 



d k 

{lo - Wa) + ka -r^ (w - 
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where the notation ka- = k{LO„) was introduced, and hence^ 

1 /"^ /,2 & , , , (Pk , ,n\ 



X Et^^ (r, w — Wct) exp(— iwt) duj + c. c. 
= {change variable of integration w' = w — Wo-} 



^i:f ( 



dk 
doj 



uj' + kc 



fk_ 



use 



X (r, w') exp(— iw't) duo' exp(— iwo-i) + c c. 



/oo 



du 



d__ d^k 



^^|^)E..(r,i) + c.c. 



As this result is inserted back into the wave equation (2), one obtains 



- ^ exp(-ia;<,i) V x V x E„^(r,t) 



+ 2ika 



dk 
du} 



d__ d^ 
.dt "d^u^^di^ 



)Ec.<,(r,i) 



+ c. c. 



a2p(NL)(r,t) 



= -Mo 



5^1 



5^ P£'^^i^)(r,i)exp(-ia;<.t)+c.c. 



: Mo J ^ p(NL) (j.^ exp(-ia;,t) + c. c. 



w„>0 



As we separate out the respective frequency components at a; = of this equation, one obtains 
the time dependent wave equation for the temporal envelope components of the electric field as 



VxVxE„„(r,i)-(A=2+i2fc,l|-fc<.0L ^)E..(r,t) = Mo^M'!^Hr,i), (3) 



where 



_/dk \-i 
^ \d^ ui^J 



Notice that unless we apply the second approximation in the Taylor expansion of k^ {lo) , terms 
containing the squares of the derivatives will appear, which will lead to wave equations that differ 
from the ones given by Butcher and Cotter. In particular, this situation will arise even if one uses 
the suggested expansion given by Eq. (7.23) in Butcher and Cotter's book, which hence should be 
taken with some care if one wish to build a strict foundation for the time-dependent wave equation. 
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Three practical approximations 

[1] The infinite plane wave approximation, 

E<^^(r,i) = E<^^(^,t)±e, 
[2] Unidirectional propagation, 



V X Vx 



Eo,,, {z, t) = A^^ {z, t) exp(±ifc<,2) 



(9^ A (9 A 

for waves propagating in the positive /negative ^;-direction. In this case, the real- valued electric 
field hence takes the form 



E(r,t)= Re[E„<,(r,t)exp(-iw<,t)] 

= ^ Re[At^_^ {z, t) exp{±ikaZ - iu>at)] 

= XI |Aa;„(2,i)|Re{exp[iA;^2:=F«w<^i + «</'(2^)]} 

= ^ \A.uj,{z,t)\cos{kaZTL^<7t + ^{z)), 



u>„>0 



where 4'{z, t) describes the spatially and temporally varying phase of the complex-valued slowly 
varying envelope fimction A^^{z,t) of the electric field. 



[3] The slowly varying envelope approximation 

dz^ 



dz 



These approximations, whenever applicable, further reduce the time dependent wave equation to 

- 50L.|i)^-(-'> - -^Pr'(M)e.p(.....). (4) 

This form of the wave equation is identical to Butcher and Cotter's Eq. (7.24), with the exception 
that here waves propagating in positive (upper signs) as well as negative (lower signs) z-direction 
are considered. 
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Monochromatic light 

Monochromatic optical field 

E(r, t) = Y, Re[E„„ (r) exp(-iw<,t)], > 

a 

E(r, a;) = ^ ^[E„^ {v)6{u - uj,) + E*^ (r)<5(u; + 

(T 

Polarization density induced by monochromatic optical field 

P(")(r, t)= J2 R'e[P£"j exp(-iw<,t)], = + uj2 + ■ ■ ■ + uJn 

(For construction of p£"^ , see notes on the Butcher and Cotter convention handed out during the 

third lecture.) 

Monochromatic light - Time independent problems 

For strictly monochromatic light, as for example the output light of continuous wave lasers, the 
temporal field envelopes are constants in time, and the wave equation (3) is reduced to 

V X V X E„„ (r) - klE^^ (r) = moo^^pL?) (r) . (5) 

By applying the above listed approximations, one immediately finds the monochromatic, time 
independent form of Eq. (3) in the infinite plane wave limit and slowly varying approximation as 

where the upper/lower sign correspond to a wave propagating in the positive/negative z-direction. 
This equation corresponds to Butcher and Cotter's Eq. (7.17). 
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Example I: Optical Kerr-efFect - Time independent case 

In this example, we consider continuous wave propagation^ in optical Kerr-media, using light 
polarized in the x-direction and propagating along the positive direction of the 2-axis, 

E(r,t) = Rc[Etj(z) cxp(— itjt)], Ejj(z) = A^{z) cxp{ikz) = exA'^{z) cxp{ikz), 

where, as previously, k = uno/c. From material handed out during the third lecture (notes on the 
Butcher and Cotter convention), the nonlinear polarization density for ar-polarized light is given 
asPr)=p(f\with 

p(J^=eo{3/A)e,X^^Li-^;u,,u,,-u,)\E;^fEZ 
= £o(3/4)xEL|Ec.|'E<, 
= eo(3/4)x£LL|At.pA„exp(iA;2;), 
and the time independent wave equation for the field envelope A^^, using Eq. (6), becomes 

= 5o(3/4)xilL|A^PA^ exp(iA:z) exp(-ifcz) 

-i^yis) lA |2a 
= {since k = Ljno{Lo)/c} 
-i — y(3) IA I^A 

= {in analogy with Butcher and Cotter Eq. (6.63), n2 = (3/8no)Xxxxx} 

C 

or, equivalently, in its scalar form 

If the medium of interest now is analyzed at an angular frequency far from any resonance, we may 
look for solutions to this equation with |A(^| being constant (for a lossless medium). For such a 
case it is straightforward to integrate the final wave equation to yield the general solution 

A^{z) = A^{zo)exp[iwn2\A^{zo)fz/c], 
or, again equivalently, in the scalar form 

^Ziz) = AZizo) oxp[tujn2\AZ{zo)\'^z/c], 

which hence gives the solution for the real-valued electric field E(r, t) as 

E(r,t) = Re[E^(z)exp(-iwt)] 

= Re{Ai^{z) exp[i(fc2: — cot)]} 

= Re{At^(^;o) exp[i(fc2: + ujn2\Ai^{zo)\'^z/c- u)t)]}. 

From this solution, one immediately finds that the wave propagates with an effective propagation 
constant 

k + ujn2\A^{zQ)f/c = {uj/c){na + n2|A^(zo)|^), 
that is to say, experiencing the intensity dependent refractive index 

"eff = "0 +"2|Aa;(2;o)|^- 



^ That is to say. a time independent problem with the temporal envelope of the electrical field 
being constant in time. 
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Example II: Optical Kerr-effect - Time dependent case 

We now consider a time dependent envelope E^^{z,t), of an optical wave propagating in the same 
medium and geometry as in the previous example, for which now 

E(r, t) = Re[Ei^{z, t) exp{—iLot)], E^^z, t) = A^^^z, t) exp{ikz) = e^A^^z, t) exp{ikz). 

The proper wave equation to apply for this case is the time dependent wave equation (4), and since 
the nonlinear polarization density of the medium still is given by the optical Kerr-effect, we obtain 



is the starting point for analysis of solitons and solitary waves in optical Kerr-media. The obtained 
equation is the non-normalized form of the in nonlinear physics (not only nonlinear optics!) often 
encountered nonlinear Schrddinger equation (or NLSE, as its common acronym yields). 

For a discussion on the transformation that cast the nonlinear Schrodinger equation into the 
normalized nonlinear Schrddinger equation, see Butcher and Cotter's book, page 240. 




= £o(3/4)xEL|A„(^, t)fA^{z, t) eMikz) exp(-ifc0) 

ZK ^ ' 




The resulting wave equation 
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